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Chiral meson masses at the critical point 
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Chiral meson masses at the critical point is investigated using QCD in the improved 
ladder approximation. We calculate the effective potential at finite temperature T and 
quark chemical potential /r and find the critical point at T ~ 95 MeV, p ~ 290 MeV. The 
chiral meson masses are determined from a second derivative of the effective potential 
at its minimum. We find that the sigma goes to massless at the critical point while 
pion remains massive. Our results are consistent with that of the linear sigma model in 
contrast to the Nambu-Jona-Lasinio model. 
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Recently there has been great interest in studying the phase structure of quantum chromodynamics 
(QCD). We expect that at sufficiently high temperature and/or density the QCD vacuum undergoes a 
phase transition into a chirally symmetric/deconfinement phase 0, a color superconducting phase ||]. 
They may be realized in high-energy heavy-ion collisions at the BNL Relativistic Heavy Ion Collider 
(RHIC) and CERN Large Hadron Collider (LHC). These phase transitions are also important in the 
physics of neutron (or quark) stars and the early universe. 

In massless two-flavor QCD, as confirmed by using the lattice simulation and several effective 
theories the chiral phase transition at high ternperature is probably second order. On the other 

hand, at high density a first order one is expected [|4|-|l0|. These observations indicate the existence of 
a tricritical point in the phase diagram for zero current quark masses and a critical point, where 
the first order transition ends, for nonzero current quark masses. The existence and the location of 
the critical point has been studied also by using the recently proposed lattice QCD method 0- It is 
expected that the sigma field becomes massless at the critical point while the pion field remains massive 
and the critical point is in the same universality class as the three dimensional Ising model [ p^ . Recently, 
it has been proposed that this point may lead to characteristic signatures which enable us to explore the 
phase structure of QCD in heavy-ion collisions [^. The behavior of the chiral meson masses around 
the critical point is studied within the framework of the linear sigma model and Nambu-Jona-Lasinio 
(NJL) model [^. In the linear sigma model, the sigma mass goes to zero at the critical point. However, 
in the NJL model with the random phase approximation the sigma remains massive at the critical point. 
The chiral meson masses at finite temperature and density has been studied in another QCD motivated 
model however, their behavior around the critical point has not been known yet. In this paper we 
investigate where the critical point locates and how the chiral meson masses behave around the critical 
point using the so-called QCD in the improved ladder approximation ra¬ 
in the following we fix the scale parameter of our model by the condition Aqcd = 1 except for 
numerical calculations. Its value is determined by the condition /,r = 93 MeV at T = ^ = 0. 

At zero temperature and zero chemical potential, the Cornwall-Jackiw-Tomboulis effective potential 

f for QCD in the improved ladder approximation is expressed as a functional of T,{pe) and Y,^{pe) 
the scalar and the pseudo-scalar part of the dynamical mass function of the quark respectively 
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is introduced, pe denotes the Euclidean momentum, g^{pE) is the QCD running coupling of one-loop 
order, C 2 is the quadratic Casimir operator for color SU (fVc) group, and an overall factor (the number 
of light flavors times the number of colors) is omitted. Note that in the derivation of Eq. (0), the 
Higashijima-Miransky approximation has been used and an infrared hnite running coupling and 

quark mass functions like Eqs. (i,® and (|^) are assumed. In this paper, we use the following effective 
running coupling [p^ 
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where pn is an infrared regularization parameter to keep the QCD running coupling from blowing up at 
Pe = ^{^qcd), and a = 8/9 (we use three-flavor, three-color effective running coupling). Correspond¬ 
ing to the above running coupling, the Schwinger-Dyson equations, which are the extremum conditions 
of Eq. (0) with respect to ^{pe) or E 5 (p£;), suggest the following trial mass functions 
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where a (cts) is related to the renormalization group invariant (RGI) quark condensate (qq) {{qi^zq)) 
as tJ = 27r^a(gg)/3 (us = 2Tr'^a{qij^q)/3) and mu is the RGI current quark mass. The explicit chiral 
symmetry breaking term is introduced in Yi(jje) alone as it should be. 

We now turn to discussions about the effective potential at hnite temperature and density. In order to 
calculate the effective potential at hnite temperature and density we apply the imaginary time formalism 
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where = (2n -|- l)7rr (n = 0, ±I,±2,---) is the fermion Matsubara frequency and p, represents 
the quark chemical potential. In addition, we need to dehne the running coupling and the trial mass 
functions at hnite temperature and density. For simplicity, we adopt the following functional form for 
2?T,/i(p), Sr,/i(p), and T, 5 .^T,fj.ip) by replacing p 4 in V(j)e), ^(pe), and S 5 (p_B)with LUn. 
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A few comments are in order. 

In Eq. (0) we do not introduce the p dependence in 'DT,ti(p)- The gluon momentum squared is the 
most natural argument of the running coupling at zero temperature and density, in the light of the 
chiral Ward-Takahashi identity [^,^. Then it is reasonable to assume that ‘Dr.fj.ip) does not depend 
on the quark chemical potential. In addition, the screening mass is not included in Eq. ( 0 ). As concerns 
the mass functions, we use the same function as Eqs. (0) and (0) except that we replace p 4 with 
The quark wave function does not suffer the renormalization in the Landau gauge for T = p = 0, while, 
the same does not hold for hnite temperature and/or density. Furthermore, we neglect the T and/or p 
dependent terms in the quark and gluon propagators that arise from the perturbative expansion. 

Substituting Eqs. (0), 0 and (|) into (0) and considering the differentiation with respect to p'Y to 
be that with respect to |p", we can write down the effective potential V{a,a^]mR) (see the Appendix 
of Ref. §). 

In numerical calculations, since it was known that the temperature and chemical potential dependence 
of the quantities such as {qq) and are stable under the change of the infrared regularization parameter 










0. Therefore, in the first place, we fix ln(p|j/AQ^^) = 0.1 and determine the value of Aqcd by the 
condition = 93 MeV at T = /r = 0 and in the chiral limit; i.e., mn = 0. In this case, the pion decay 
constant is approximately given by Pagels-Stoker formula : 
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and we obtain Aqcd = 738 MeV. Secondly, we assume the light quarks {u and d) are degenerate in mass 
and take the current quark mass evaluated at the renormalization point k = 1 GeV as m„(lGeV) = 
7MeV. Using the one-loop evolution formula, the RGI current quark mass rriE extracted from the 
above-mentioned value becomes to_r = 7.6 x 10 “^Aqcd. A discussion of the effective potential and the 
chiral symmetry restoration at high temperature and/or density can be found in Ref. |^. 

Figure 1 shows the phase diagram in the chiral limit in ns-T plane. The baryon number density ns 
is defined as ub = nq/3 = —{l/3)dV/dp, with V renormalized so that it has the correct free theory 
behavior at cr = <75 = 0. At T = 0, for example, there is a mixed phase that consists of massive quarks 
with ^ = 1.5no and massless quarks with = 4.3no where uq = 0.17fm~^ is normal nuclear 
matter density. The phase diagram for 77 iii(lGeV) = 7MeV case is shown in Fig. 2. The critical point 
E where the first order phase transition ends is found at pE — 290 MeV, Te — 95 MeV; 


~ 0.64, cx 0.70, (11) 
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where Pcrit is the critical chemical potential at T = 0 , and Tcross is the temperature where is 
minimized and starts to increase at /r = 0. We have confirmed that by the finite current quark 
mass the critical point is moved from the tricritical point, which we have found in the previous paper 
1^ , toward larger value of p and smaller value of T. 

The values of the chiral meson masses are defined by multiplying the second derivative by the appro¬ 
priate factor /: 
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Here, “min” at the end of the equations means that they are evaluated at the minimum of V(a, us; to_r). 
In this paper, we do not examine the factor /; rather, we fix Mt^{T = p = 0) to 140 MeV, then, 
Mcr{T = p = 0) turns out to be 668 MeV. The p dependence of and at T = Tg is shown in 
Fig. 3. We find that sigma almost goes to massless at the critical point, while pion remains massive. 
The chiral meson masses at p = pe show behavior similar to Fig. 3 as functions of T. These behavior 
is consistent with that obtained in the linear the sigma model Q . 

In conclusions, we investigated the phase structure and the chiral meson masses at the critical point 
from QCD in the improved ladder approximation. Using the variational approach, we calculated the 
effective potential at finite temperature and chemical potential. In the chiral limit, we found the mixed 
phase that consists of a low density chirally broken phase and a high density chirally symmetric phase 
when T is sufficiently small. For explicit chiral symmetry breaking case, the critical point was found at 
Pe — 290 MeV, Te — 95 MeV. The value oi pe seems to be too large in light of the experiments at 
RHIC and LHC. By including strange quark mass, however, it would be reduced. The meson masses 
was obtained from the curvature at the minimum of the effective potential. We found that M„ almost 
goes to zero at the critical point while Mjr remains nonzero. The result is consistent with that obtained 
within the framework of the linear sigma model in the mean-field approximation [^. We have also 
confirmed that {qq) at the tricritical point behaves as 
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Our results can be understood by the Landau-Ginzburg effective potential and partly support the 
experimental signature proposed in Ref. |13|. 

Finally, some comments are in order. In this paper we did not take into account the screening of 
the gluon, that is to say, the Debye screening for the electric gluons. Furthermore, it has been shown 
that at finite temperature the wave function renormalization makes non-negligible contribution to the 
critical line and the behavior of the order parameter 0 ]- They may affect the precise location of the 











critical point although the main feature of this work might not change. In any case, it is preferable to 
take into account the screening of the gluon, wave function renormalization and the effects of s quark 
in the future work. 
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FIG. 1. The phase diagram in the chiral limit as a function of the baryon density and the temperature. The 
solid line indicates the phase transition of second order and the dashed lines indicate that of hrst order. 
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FIG. 2. The phase diagram in fi-T plane. The solid line indicates the phase transition of second order and 
the dashed line indicates that of first order. The thin line corresponds to the chiral limit and the thick line 
corresponds to mi{(lGeV) = 7MeV case. The points P and E represent the tricritical point and the critical 
point, respectively. We note that the critical chemical potential at T = 0 is slightly (about 2%) larger than 
previous papers [9] because we have improved the numerical calculation. 










